HUPD-9511 
March, 1995 



Fermionic excitation in quark matter 



Kazuhiro Tanaka 

Department of Physics, Hiroshima University, Higashihiroshima-shi, Hiroshima 724, Japar^ 



Abstract: We discuss the fermionic excitation spectra of the dense QCD and QED 
plasma at zero temperature, and analyze the stability of the normal (perturbative) ground 
state. The standard re-summed perturbation theory breaks down because the infrared (IR) 
divergences show up in the spectra near the Fermi surface. We employ the effective field 
theory approach and the renormalization group to investigate the dynamics near the Fermi 
surface. Our results indicate that the dense QCD and QED plasma are non-Fermi-liquids. 



Talk presented at YITP Workshop on 

"From Hadronic Matter to Quark Matter: Evolving View of Hadronic Matter' 1 . 
YITP, Kyoto Japan, October (1994). 

(to be published in Supplement of Progress of Theoretical Physics) 



E-mail: tanakak@theo.phys.sci.hiroshima-u.ac.jp 



1. Introduction 

There are several situations where hot and/or dense deconfined quark-gluon plasma (QCD 
plasma) might be formed; e.g., the interior of the neutron stars, the early universe, and the 
relativistic heavy-ion collisions. It is commonly accepted that for the QCD plasma the free 
gas description gives a useful starting approximation. The another example for which such 
treatments are frequently used is the QED plasma (electron gas). This corresponds to the 
following assumption: for sufficiently hot and/or dense case, perturbation theory (PT) along 
with some appropriate resummation procedure (re-summed perturbation theory, rsPT) gives 
a successful scheme to describe the dynamics. 

In this talk we concentrate on the dense case, by which we mean the high-density but 
zero-temperature case. For this case, the above assumption corresponds to that the ground 
state is the normal Fermi-liquid state. We discuss the fermionic excitation spectra in the 
dense plasma, and examine the validity of the assumption. We emphasize that the physics 
of the dense QCD and QED plasma is profoundly connected with the IR behavior of the 
theory. 

2. Perturbation theory calculation of the ground-state energy 

A peculiar property of QCD is the asymptotic freedom, which is a basis of successful 
applications of perturbative methods to high-energy processes!]]. Naively, one would expect 
that the asymptotic freedom would support PT starting from the free quark gas. Such a 
picture was pushed forward by Collins and Perry 0. They employed a renormalization group 
(RG) argument similar to those in the high-energy process. The following behavior of the 
thermodynamic quantities, when the Fermi momentum pp is scaled up by k, was obtained: 

X(Kp F ,g,m,n) = k Dx X(p F ,g(K),m(K),(i), (1) 

where X stands for some thermodynamic quantity like the energy density, and Dx is its mass 
dimension, g and m are the gauge coupling constant and the quark mass renormalized at the 
scale /i, while g(n) and m(/«) are the corresponding running quantities. Because g~(K) — > 
as k — ► oo, (TTJ) seems to support the naive expectation discussed above. 

There is a subtle point in this sort of RG arguments: the difficulties might arise from 
the IR divergences, because the gluons are massless and also m(/c) —>■ as k — > oo. In the 
medium, however, the gluons acquire the Debye screening mass and the plasma mass due to 
the interaction with the quark gas||]. Collins and Perry argued that these "effective masses" 
would act as an IR cut-off. 

Following this line, and in order to demonstrate explicitly IR safety of PT at high density, 
the calculation of the ground-state energy, up to and including the effects of 0(a 2 ) (a = 
g 2 /4:7i), was carried out in (|, At 0(a 2 ), the diagrams, which involve the 3-, 4-point gluon 
couplings or the ghost-gluon coupling, contribute. As for the IR divergences, their results can 
be summarized as follows:^ the second order (0(a)) graph is free of the IR divergences. At 
the fourth order 0(a 2 ), the individual diagrams contain the IR divergences, but most of them 
cancel out between the different diagrams. The only one exception is the diagram involving 
the two fermion bubbles; thus the order by order PT is suffered from the IR divergences. 
Fortunately, however, the resummation of the higher order graphs of the RPA-type (ring 

2 The ultraviolet (UV) divergences can be completely removed by the counter terms determined at pp = 0. 
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diagrams) renders the result finite. This is similar to the well-known nonrelativistic (NR) 
electron gas calculation ||. The resummation of the RPA produces a finite static screening 
mass to the "dangerous modes" of the gluons, and the difficulty can be avoided. The rsPT 
for the dense QCD plasma is IR safe in accord with the conjecture of ref. @, at least for the 
computation of the ground-state energy. 

There is one price due to the resummation: The re-summed RPA series (the sum of the 
ring diagrams involving two or more fermion bubbles) give the contribution of 0(a 2 loga), 
which are dominant compared with the other 0(a 2 ) graphs. 

3. Fermionic excitation spectra 

In this section we compute the quark excitation spectra in the dense QCD plasma in the 
framework of rsPT discussed in sect. 2. 

The excitation energy of the quarks in the medium can be defined as the change in the 
total energy when a quark is added (removed) into (from) the system. This can be expressed 
as a functional derivative of the energy density S with respect to the Fermi distribution 
function n(p): 

<p) = /fr, (2) 

e(p) is nothing but the quasiparticle energy in the Landau Fermi-liquid theory 0. In PT, 
the (free) quark propagator in the rest frame of the T = plasma is|3| 

s{ p) = + ^-^rr 5 (p° ~ e(p)Hp) = s f (p) + s D ( P ), (3) 

p — m + is E(p) 



with E(p) = \Jp 2 + m 2 and n(p) = 6(p F — |p|). Because Sp involves n(p), the functional 
derivative (|2p cuts the one quark-line in the diagrams for the energy density. Including the 
contributions up to 0(a 2 log a) for S, we obtain as the corrections to the free quark spectrum 
E(p): 

Tfl 

ei(p) = -j^u{p)E(p)u(p). (4) 

u(p) is the free quark spinor normalized asuu= 1. S(p) is the RPA-type self-energy given 
by the diagrams of fig.l. For the computation of e\(p), we perform the Wick rotation for the 
gluon loop momentum integration. The procedure is similar to the one discussed in 0, ||. 
One obtains the two types of contributions, depending on whether the intermediate quark 
line connecting the external legs in the diagrams of fig.l is off-shell or on-shell. We denote 
the off-shell part by eY(.P) while the on-shell part by ef (p), as £\{p) = e^ip) + € i\p)- e i can 
be compactly expressed as 

ef (p) = 2N,N e £' |^ £dcos0/(p )9 ). (5) 

cos8 = p ■ q/pq, and f(p, q) is the exchange-type interaction (averaged over spin and color) 
between the two quarks with the momenta p, q. Here and in the following, we consider the 
case of SU(iV c ) color and the Aj-light quarks. 

eY is given as a (Wick-rotated) 4-dimensional euclidean integral. It can be shown, by 



explicit computation, that is free of the IR divergences |10]. On the other hand, e*p is IR 
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divergent! / assumes the following form| 

f<P,Q)= E 



A 



x 



X=L,T 



(p - q) 2 - n x (p - q)' 



(6) 



where "L" and "T" denote the longitudinal and the transverse modes of the gluons. They 
are the two independent gluonic modes in the medium. IIx is the self-energy for these 
modes due to the RPA resummation. It is given by the gluon self-energy 11^ as IJi,{k) = 
nu(k) — II 00 (k), IlT{k) = n^{k) if we choose k to lie on the x-axis. Ax is the dimensionless 
function accounting for the vertex structure of the quark-gluon coupling and is proportional 
to a. The IR divergence shows up near the Fermi surface. For p = pp, we obtain 



ef(p) (p-p F ) 



NfNc 
4vr 2 



E 

X=L,T 



1 A x 
dcosO - . 

-i 1 - cos 9 + n x (p - q) 



(7) 



and n x {k) = 
™™ a >0. 



where the integrand is evaluated at the Fermi surface = \q\ = pp) 
(n x (k)/2p F ) k0=(} . The integral would be logarithmically divergent if IIx = 0{9 2 ) as 
This is the IR divergence due to the soft gluons with the momentum p M — g M — > 0. 

Now the key quantity which determine the IR behavior of the fermionic excitation spectra 
is the so-called "screening mass" of the gluons defined by 

M\ = lim lim n x (k, k ). (8) 

k—>0 ko—>0 

In the RPA resummation scheme, IIx contains the quark one-loop diagrams. One finds: 



a 



Mi = 2N f -pU 

7T 



Ml 



0. 



(9) 



These are the gauge invariant results. Therefore, the L-modes give the IR safe result, while 
the T-modes cause the IR divergences. The screening of the L-modes corresponds to the 
Coulomb screening in the NR electron gas||. The non-screening of the T-modes (magnetic 
modes) is a result of gauge invariance k fJi II fll/ = in the framework of PT. 
We now find for the quark excitation spectra: 



e(p) = £f + (p - Pf)v + 0[{p- p F ) 



where ejr is the Fermi energy, and 



V = Vq - C F V n 



7T 



log 



Pf 



P~Pf 



0(1) 



(10) 



;n) 



with t> = Pf/E(pf),Cf = (N 2 — 1)/2N C . «(pf) is the running coupling constant at the 
scale pp. The result indicates the logarithmic singularity at p = pf, and thus the negatively 
infinite slope of the spectra at the Fermi surface (vf = de{p)/dp\ p=PF — > — oo). The behavior 
of the spectra is shown in fig. 2. Several comments are in order here: 

(a) Also for the case of the hot QCD plasma (at pp = 0), the T-modes of the gluons do not 
get screened, at least up to 0(a)||: it was pointed outJTl] that this causes the breakdown 
of PT for the thermodynamic potential at the sixth order. On the other hand, for the case 
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of the dense plasma, there has been no clear indication of such difficulties: all energy density 
computations by rsPT have been worked out successfully as discussed in sect. 2. Our results 
reveal that the RPA breaks down if one computes the excitation spectra. 

(b) In the present treatment, the non-abelian character of the theory does not enter 
explicitly, except for the running coupling constant a{pp) and the color factors of flTTI). 
Thus the result can be easily translated into the case of an abelian gauge theory, i.e., into 
the QED plasma. This indicates that a naive extension of the famous work on the NR 
electron gas by Gell-Mann and Brueckner|6|] to the relativistic case breaks down:[] in the 
NR case, the quasiparticle energy in the RPA scheme was obtained successfully by Gell- 
Mann[TJ|]. Our results (|H]),([11|) do not reduce to Gell-Mann's result even if we go to the NR 
limit {ppjra — > 0). The discrepancy comes from the IR divergent contributions discussed 



above. Our results confirm the argument of [|14]] by employing the full relativistic framework. 



(c) The breakdown of rsPT poses a question on the usefulness of the free gas description 
of the dense QCD and QED plasma. The spectra with the inversion may indicate the 
instability of the normal (Fermi-liquid) state (see also |TJ], [15], [16], Q). 



The breakdown of the RPA as rsPT for high density has a great impact. Before drawing 
physical implications, however, we have to investigate the effects beyond the RPA, and 
examine whether those effects cancel the above IR divergences or not. In the next section 
we employ the RG approach and try to go beyond the RPA. 

4. Renormalization group approach 

In this section we investigate the ground state of the dense QCD plasma in the framework 
of the renormalization group (RG) a la Wilson. The Wilson RG is a useful tool to analyze the 
low energy dynamics of the system and to construct an effective field theory. It is especially 
powerful when the low energy dynamics is dominated by gapless (massless) excitations. In 
this case, the effective field theory corresponds to an IR fixed point of the RG transformation. 



A well-known example is the application to critical phenomena] 17 



Recently, there has been progress in the applications of the RG to (NR) many-fermion 



systems|18], [Hj, gJJ. The key observation is that the system of many-fermions in the 
normal state is dominated by gapless excitations at low energy, due to particles and holes 
near the Fermi surface. The Landau Fermi-liquid theory corresponds to the Gaussian (mean- 
field) fixed point. This view point also opens a possibility to classify the possible phases of 
the system by a few number of "relevant perturbations" which would lead to a RG flow 
toward another IR stable fixed point rather than the Fermi-liquid fixed point. For example, 
it has been shown Hl8j |19|, |20| that, for the case of 3-dimensional NR Fermi systems with 
rotational invariance, only the pairing interactions leading to the Cooper pairs could be 
(marginally) relevant; otherwise, the system is described as a normal Fermi-liquid. 

Now we apply the RG approach to the dense QCD plasma, and investigate whether the 
Fermi-liquid fixed point is a IR stable one or not. If it were not stable, the assumption 
discussed in sect.l would be wrong. As the first step, we have to write down the low energy 
effective action S\. It is obtained by integrating out the "fast modes", which have the 
momenta greater than the UV cut-off A, in the original action. The modes within A, which 
we call the "slow modes" , constitute our "physical" degrees of freedom dominating the low 



3 The ground-state energy of the relativistic electron gas was obtained successfully in|E 
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energy dynamics, and S\ should be expressed by them. One important point is that all 
possible relevant or marginal interactions (perturbations) between the slow modes, except 
those forbidden by symmetry, should be included in S\. for our RG analysis is: 

Sa = I d 4 pii(p) 1 \ P -e F -v F (p- PF )]iPip)+ [ d 4 qq 2 (A«(q)) 2 

J\p-p F \<A J J\q\<A V 7 

+ / dV / d 4 p j d 4 q {2ixf5 {i) (p'-p- q) g^{p')-i i t a ij(p)A a i (q) 

J\p'-p F \<A J\p-p F \<A J\q\<A 

+ (4-fermion interaction terms) + (gluon interaction terms). (12) 

Here the quark (ip) and the gluon (Af) fields are in the momentum space.0 The momentum 
integrations are restricted within A and involve the slow modes only. Note that, for the 
fermions, the slow modes are those near the Fermi surface. Their excitation spectra are 
linearized as e(p) = e F + v F {p — p F ). Af stand for the T-modes of the gluons. The L-modes 
of the gluons do not appear in (0): as we saw in sect. 3, the "mass" of the gluons relevant for 
the fermionic excitation spectra is the screening mass. The L-modes become massive due to 
the Debye screening, and thus decouple from the low energy dynamics. On the other hand, 
the screening mass for the T-modes is forbidden by gauge invariance. The integration of the 
L-modes and other fast modes generate the 4-fermion interactions, but only a few types of 
them, which satisfies special kinematical conditions, could be marginal (at the tree level). 
The "gluon interaction terms" may involve the 3- and 4-point couplings between the gluons, 
and also the gauge terms (the gauge fixing terms and the Faddeev-Popov ghost terms), but 
the explicit forms of them are irrelevant for the following discussion. 

As a result of mode elimination, the coupling parameters (v F , g, the 4-fermion couplings, 
etc.) depend on A, and thus are the running couplings in the RG. The partition function is 
given as the functional integral by the slow modes 

dipdtfdA] e iSA . (13) 

We perform the RG transformation and derive the RG equations for the running couplings. 
We divide our physical space further into the slow modes (\p — p F \,q < A/s, s > 1) and 
the fast modes (A/s < \p — p F \,q < A), and integrate the fast modes out in dl3|). The 
scale transformation \p — p F \,q — > \p — p F \/s, q/s recovers the phase space reduced by mode 
elimination. The comparison of the new effective action with ( |T2| ) gives a flow of the coupling 
parameters when A is scaled down. In order to trace the RG flow toward the Fermi-liquid 
fixed point, we perform the scale transformation to the quark and the gluon fields as ip —> 
s 3 ^ 2 y/Z^tp; A® — > s z \fZ~?,A\i where Z2 and Z% are the renormalization constants for these 
fields. 

The Fermi velocity v F determines the behavior of the fermionic excitation spectra near 
the Fermi surface. This is a marginal coupling at the tree level. By performing mode 
elimination at the 1-loop level (see fig.3), we obtain the RG equation: 

. dv F a 9 . . 

A-f = C F -v 2 F , 14 
dA 71 

4 These fields stand for the fermionic and the bosonic elementary excitations with the same quantum 
numbers as the quark and the gluon; we call them simply by the quark and the gluon in this section. 
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with a = g 2 /4:7T. This should be combined with the /3-function, (3 = A(da/dA). As a 
boundary condition, we match our effective field theory with the QCD for the zero density, 
at the initial (high-energy) scale A (~ Pf) of the successive RG transformation toward 
the Fermi surface: a (A) = 1/ [/? log(A/A QCD )], or /3 = -(3 a(A) + 0(a 2 ), for A > A . 
(Po = (11 N c — 2Nf)/6ir and Aqcd is the QCD scale parameter.) This implies, by the 
continuity, a (A) grows at the beginning if we scale down A from Ao- When we go to the 
lower energy, we have the two possibilities, depending on the behavior of (3(a) for a > a(A ): 

(i) (3{a) < (i.e., no zero of /3(a)) for a > a(A ). 

(ii) (3(a*) = for some a* > a(A ). 

The case (i) corresponds to the case where a is a marginally relevant coupling and 
continues to grow at low energy. In this case, the low energy effective theory is a strong 
coupling theory and thus does not correspond to the Fermi-liquid fixed point. 

The case (ii) corresponds to the case where a* gives an IR fixed point. If a* is small, 
this may correspond to the Fermi-liquid fixed point, and thus needs more detailed analysis. 
In this case, (0) can be easily integrated to give 

(/K] = M A o) n ^ 

F{ ] l + C F («*AOMAo)log(A /A)- 1 ) 

Here vp(A ) > 0. The result implies Vp — > when one goes to the low energy limit A — > 0, 
i.e., near the Fermi surface. Note that if we expand (|15|) perturbatively and retain the 
leading term, we recover the previous result (|TT|) . This indicates that the discussion of 
sect. 3 corresponds to assuming the case (ii); the RG makes it possible to sum up the leading 
logarithms in all orders of PT. 

One might say that the result is "safe" if the case (ii) is realized. We have no IR difficulty 
(|15D indicates that Vp is a marginally irrelevant coupling at low energy. So, the low energy 
effective theory would be given by the Fermi-liquid fixed point with a = a* and Vp = 0. 
However, we argue that this fixed point will never be stable: first, because a* ^ 0, there 
exist the long-range attractive interaction between the quarks due to the exchange of the 
T-modes of the gluons. Second, the density of states at the Fermi surface N(0) oc 1/vp is 
infinitely large (see fig.2)|] The infinite density of states will enhance the attractive channels 
of the 4-fermion interactions, and will drive the instability of the Fermi-liquid fixed point. 
For example, the BCS-type paring interaction V gets the 1-loop correction ~ — N(0)V 2 , and 
is strongly enhanced at low energy for V < 0. As an origin of the attractive V, we need not 
other degrees of freedom like the phonon in the usual superconductors; it could be generated 



solely through mode elimination in gauge theory, e.g., by the Kohn-Luttinger effect p2 . 

Our RG analysis indicates the two possibilities for the effective field theory for the dense 
QCD plasma: (i) a strong coupling theory with the large a; (ii) a "would-be fixed point" 
with a = a*, Vp = 0. We conclude that the ground state of dense QCD plasma is not the 
normal Fermi-liquid state for both cases. 

In order to determine which case is realized, we have to compute the /3-function from 
the relevant 1-loop diagrams. For the simpler case of the dense QED plasma, we can show 



it corresponds to (ii) at 1-loop HlQ[. (Note, in this case, (3 < 0.) Thus, the results similar 



5 In the NR case where the T-modes of the gauge fields are not included, the running of vf is not so 
significant, and vp is considered to go to a nonzero value for A — > 0fll8|]. 
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to fljjD can be obtained: the ground state of the electron gas may be the superconducting 
state, even without the phonon degrees of freedom.^ 

The behavior u.f(A) — > (A — > 0) is similar to that of a marginal Fermi-liquid.Q The 
marginal Fermi-liquid is a successful phenomenological model proposed to explain the puz- 
zling normal state properties of Cu-0 high-temperature superconductors ||25|| . That a proper 
relativistic treatment of the electron gas might lead to superconductivity of a marginal Fermi- 
liquid would be an interesting possibility as a model of high-temperature superconductivity. 

5. Conclusion 

We discussed the fermionic excitation spectra of the dense QCD (and QED) plasma. It 
is shown that the standard rsPT breaks down because the IR divergences show up in the 
spectra near the Fermi surface. The divergence is due to the non-screening of the transverse 
modes of the gauge fields. By employing the effective field theory approach and the RG, 
we classified the two possibilities for the ground state of the QCD plasma: (i) a strong 
coupling theory where the nonperturbative effects will be important; (ii) a would-be IR 
fixed point with a = a* and vf = 0, which would be unstable against the transition to 
the superconducting state. We concluded that the T = QCD (and QED) plasma is a 
non- Fermi-liquid. 

We stress that the physics of the dense QCD and QED plasma is not so simple on the 
contrary to the usual picture. Our results indicate the importance of the "soft degrees of 
freedom"; they could appear in the higher order terms in PT as well as in the nonpertur- 
bative effects. More efforts are needed to reveal roles of the soft degrees of freedom and to 
understand the dynamics of the QCD and QED plasma. 
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Figure Captions 



Fig. 1 The self-energy diagrams for the quarks in the RPA (full lines: quarks; wavy lines: 
gluons) . 

Fig. 2 The fermionic excitation spectra (dotted line: Vp > 0; dashed line: Vp = — oo; full 
line: Vp = 0). 

Fig. 3 The 1-loop diagrams which give the RG equation (fl^). The first diagram is due to the 
transverse gluon exchange, while the second one is due to the 4-fermion interactions. 
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